Brain Structure and Function

Evidence from in-silico, ex vivo and invasive animal studies suggested the function of neural tissue to be inextricably linked to its structural and microstructural organization. However, this link has never been assessed non-invasively in humans. Here, we established a simple quantitative framework linking cerebral oxygen metabolism to gray matter microstructural features in healthy humans using advanced non-invasive MRI methods. The cerebral metabolic rate of oxygen consumption (CMRO2) was quantified with calibrated fMRI, whereas microstructural properties were estimated from diffusion MRI data fitted to the Soma and Neurite Density Imaging (SANDI) model, yielding soma radius (R soma) and soma volume fraction (fsoma) metrics, from which we derived maps of soma surface area density (SAD) and soma numerical density (ND) under the assumption of a spherical soma. CMRO2 was positively associated with SAD (r = 0.29, p < 0.01), fsoma (r = 0.28, p < 0.05) and ND (r = 0.25, p < 0.05), and negatively associated with Rsoma (r = -0.23, p < 0.05) across cortical brain regions. These relationships are consistent with predictions from experimental and theoretical studies, providing the first quantitative evidence that cerebral metabolic demand is coupled to microstructure in humans. The positive association of energy consumption with SAD is consistent with the notion that both the number of synapses (energy-demanding structures), and the energy required to maintain a membrane potential, scale with soma surface area. Our findings establish a non-invasive MRI framework for investigating the microstructural determinants of brain energy metabolism in health and potentially in disease.

The choroid plexus (ChP) is a highly vascularized structure located in the brain ventricles that plays important roles in multiple central nervous system functions. ChP enlargement has been implicated in many neuropsychiatric conditions such as multiple sclerosis, depression, schizophrenia, and Alzheimer's disease. However, limitation in the accuracy of the current segmentation tools highlight the need for more reliable methods to quantify ChP volume in healthy and clinical human populations. The shape and size of the ChP vary drastically from person to person, and currently there are no large ChP datasets which allow comprehensive training of segmentation models. Here we introduce DataPlex, a collection of 700 structural MRI brain images where the lateral ventricle ChP has been manually labeled and checked by expert human raters. DataPlex is balanced and representative across a wide range of scanning parameters, demographics, and clinical profiles, and subject-level information will be available for each labeled image. DataPlex will be open source to support the creation of reliable segmentation tools for the human lateral ventricle ChP.

Abstract The importance of hypothalamic kisspeptin (Kiss1) neurons in controlling the reproductive axis is well established, while a population of extra-hypothalamic Kiss1 neurons is found in the medial amygdala (MeA) of both males and females in rodents and is less well understood. The MeA is connected with the olfactory bulbs, the hypothalamus and cortical areas. MeA Kiss1 neurons have been implicated in sexual behaviour and modulation of the reproductive axis; however, little is known about their action potential firing properties, which would influence kisspeptin release in target brain regions, nor whether their physiological properties are sex- or estrous cycle-dependent. In this study, whole-cell patch-clamp recordings were made from Kiss1 neurons in brain slices containing the MeA of adult female (estrous and diestrous) and male mice. MeA Kiss1 (Kiss1MeA) neurons were sexually dimorphic, with greater numbers in male compared to female mice and a slower change in excitability in estrous females compared with males. However, combining a range of electrophysiological parameters using principal component analysis did not reveal a distinct sex- or estrous cycle-linked phenotype, suggesting that the properties of Kiss1MeA neurons- and therefore kisspeptin release- are broadly similar. Interestingly, Kiss1MeA neurons had unique electrophysiological properties compared with anteroventral periventricular (AVPV) hypothalamic Kiss1 (Kiss1AVPV) neurons in female mice; with their reduced excitability, achieving action potential firing frequencies greater than 10Hz required a larger current input, indicating that Kiss1MeA neurons require stronger synaptic input to achieve firing frequencies that drive kisspeptin release. These findings provide new information about the physiological properties of Kiss1MeA neurons, showing subtle sex differences and profound regional differences between the MeA and AVPV.

Syntax processing impairment is a defining feature of primary progressive aphasia (PPA), particularly in the nonfluent/agrammatic variant. However, current models are derived largely from Indo-European languages, and it remains unclear whether the cognitive and neuroanatomical mechanisms of syntax processing generalize to Chinese, an analytic, pro-drop language that encodes grammatical relationships primarily through word order and function words rather than inflectional morphology. Participants were recruited from the international Chinese Language Assessment in PPA cohort and included 64 individuals with PPA (19 logopenic, 15 nonfluent/agrammatic, 15 semantic, and 15 unclassifiable/mixed) and 60 cognitively healthy controls. Participants completed sentence comprehension and production tasks with systematically manipulated syntactic features. Item-level linguistic determinants were examined using ridge regression, diagnostic performance using receiver operating characteristic analyses, and neuroanatomical correlates using voxel-based morphometry of structural MRI. Performance on sentence comprehension and production was influenced by grammatical and lexical features. Subject-relative clauses, passive constructions, a higher proportion of function words, longer character length, and greater lexical content were associated with poorer performance among participants with PPA. Sentence comprehension was impaired in the logopenic, nonfluent/agrammatic, and semantic variants, even when isolated to the syntactic component of the sentence comprehension task. Similarly, sentence production was impaired across all PPA variants (all P < 0.01 versus controls). No significant differences were observed among PPA variants. Adjusted receiver operating characteristic analyses demonstrated high discrimination between each PPA variant and controls, with areas under the curve ranging from 0.87 to 0.97 for sentence comprehension-syntax and from 0.91 to 0.97 for sentence production. However, pairwise comparisons showed no significant differences in discriminative performance among PPA variants after correction for multiple comparisons. Both sentence comprehension-syntax and sentence production showed associations with semantic and repetition measures in addition to motor speech measures. Voxel-based morphometry showed that higher sentence comprehension-syntax scores were associated with greater grey matter volume in left frontotemporal and cingulate regions, whereas higher sentence production scores were additionally associated with greater gray matter volume in frontostriatal structures, including the caudate, putamen, and thalamus. Syntactic processing impairment in Chinese-speaking individuals with PPA extends beyond the nonfluent/agrammatic variant and is associated with distributed frontotemporal and frontostriatal regions, as well as lexical-semantic and auditory verbal short term memory mechanisms. These findings suggest that the cognitive and neural manifestations of syntactic impairment may be shaped by language-specific typological properties, underscoring the need for linguistically informed approaches to the assessment and interpretation of PPA.

Neuroimaging research focusing on subcortical brain regions showed that they play a role in a wide range of cognitive functions and associated their structural alterations to various psychiatric disorders. Methods to analyse such structures were developed to observe changes, not simply in terms of volumes, but in terms of surfaces, which are able to detect subtle and local changes in shape and geometrical configuration. Such advanced analyses have however been largely restricted to command line environments with relatively limited access to Python users or researchers with low technical expertise. We present SubCortexMesh as a user-friendly toolbox which covers automated surface estimation from popular subcortical volume segmentations (FreeSurfer and the Functional Magnetic Resonance Imaging of the Brain Software Library (FSL)), computes shape-related vertex-wise metrics (thickness, surface area and curvature) for whole cohorts and includes statistical analyses. This is all done inside Python, automatically computing all subjects of a given preprocessed directory, with the explicit intent to minimise steps and manual coding required from the users. SubCortexMesh includes statistical tools which allow conventional random field theory-based cluster analyses on native subject metrics, standardised within common surface templates, thus containing a whole workflow necessary to run an up-to-date shape-wise subcortical analyses in one package.

Axonal projections and precise synapse formation are essential for the establishment of functional neural circuits, and their dysregulation contributes to various psychiatric and neurodegenerative disorders. Although numerous extracellular cues secreted from surrounding tissues regulate axon guidance and synapse formation, the mechanisms by which neurons locally regulate terminal patterning within their projection fields remain poorly understood. Bone morphogenetic protein 4 (BMP4) is a classical morphogen that is involved in early neural tube patterning during development. However, its postnatal function remains unclear. Previous in vitro studies using hippocampal glutamatergic neurons demonstrated that BMP4 can act locally on axons to suppress synapse formation, raising the possibility that BMP4 regulates terminal differentiation after target innervation. Here, we investigated the role of BMP4 in postnatal cholinergic circuit development using cholinergic neuron-specific Bmp4 conditional knockout (Bmp4 cKO) mice, with adult analyses performed in males. Basal forebrain cholinergic neurons project extensively throughout the neocortex and are essential for higher brain functions, including arousal and attention. Bmp4 cKO mice exhibited increased cholinergic axon terminals in the neocortex, hippocampus, and amygdala, as well as increased axonal ramification in the neocortex without affecting cytoarchitectonic organization or cholinergic neuron numbers in the basal forebrain. These findings demonstrate that BMP4 functions as a negative regulator of cholinergic terminal formation and axonal arborization, suggesting that BMP4 may act locally to restrict cholinergic axonal expansion during neural circuit maturation.

Dysfunctions in reward circuits, including projections from the orbitofrontal cortex (OFC) to the nucleus accumbens (NAcc), have been implicated in psychopathology. The OFC is functionally divided into the medial OFC (mOFC), which processes rewards, and the lateral OFC (lOFC), which responds to non-reward or punishment stimuli. We investigated the structural connectivity of the lOFC-NAcc and mOFC-NAcc pathways and their associations with psychopathological traits using tractography-defined ROI analysis. Diffusion-weighted imaging data from 165 participants (aged 22 - 35, 53% male) in the Human Connectome Project (HCP) were analysed using spherical deconvolution tractography. MegaTrack, a novel tool, enabled tract segmentation via a single manual dissection mapped back to native space. Virtual dissection identified four bilateral pathways connecting the lOFC and mOFC to the NAcc, from which fractional anisotropy (FA) was extracted, and an lOFC/mOFC ratio was computed. Psychopathological traits were measured using HCP Achenbach Adult Self-Report scores. The lOFC-NAcc pathway showed positive associations with psychopathology and the mOFC-NAcc pathway negative associations. The lOFC/mOFC ratio emerged as a stronger predictor than individual pathways, with significant associations for depression, somatic problems, anti-social personality traits, and inattention. Effects were lateralized to the left hemisphere. These findings show that the lOFC-NAcc and mOFC-NAcc pathways can be structurally separated with tractography and display opposing correlations with psychopathological traits in a healthy population. Their ratio may index reward-processing dysfunction, and future research should test whether it is sensitive to dysfunction in clinical populations.

Identifying regional brain atrophy is critical for understanding neurodegenerative disorders. Traditionally, group comparisons rely on frequentist statistics, utilising $p$-values and point estimates that offer limited characterization of uncertainty. Bayesian methods provide an alternative by estimating effect sizes and credible intervals while incorporating prior knowledge and uncertainty. This study evaluated the feasibility of Bayesian hierarchical modelling to assess regional brain atrophy in Progressive Supranuclear Palsy (PSP) and Corticobasal Syndrome (CBS) compared to healthy controls. Participants meeting diagnostic criteria for CBS or PSP, together with healthy controls, were recruited from the Cambridge Centre for Parkinson-plus. Grey matter volume, cortical thickness, and subcortical volumes were derived using FreeSurfer for 210 cortical and 36 subcortical regions defined by the Human Brainnetome Atlas. Participants with PSP were classified into Richardson syndrome (PSP-RS), corticobasal syndrome variant (PSP-CBS), frontal variant (PSP-F), and other presentations (PSP-Other). Bayesian hierarchical ANOVA models with appropriate priors were used to estimate group differences. Regions of practical equivalence (ROPE) were calculated to assess evidence of meaningful group effects. A total of 428 people were included (113 PSP-RS, 30 PSP-CBS, 23 PSP-Other, 19 PSP-F, 81 CBS, and 162 healthy controls). For grey matter volume, there was evidence of whole brain volume loss compared to controls for all disease groups except PSP-Other, which showed evidence of no difference. For cortical thickness, there was evidence of global cortical thinning in CBS compared to controls, but for no other group. For subcortical volume all groups showed evidence of volume loss compared to controls. Patterns of regional atrophy were identified in keeping with previous literature, with the largest effect observed in the left pre-motor thalamus (30\% reduction in CBS). There was greater loss of subcortical volume in PSP-RS compared to other PSP subgroups. Bayesian hierarchical modelling provides a feasible framework for quantifying regional brain atrophy, with interpretable effect size estimates and uncertainty measures. This approach complements conventional statistical methods, particularly in the assessment of rare neurodegenerative and neuropsychiatric disorders.

The mesolimbic dopamine system supports broad aspects of adaptive behavior including salience, learning, memory, and decision making, yet the structural organization supporting these functions remains incompletely mapped in humans. Using high-resolution 7T diffusion magnetic resonance imaging (up to N=173), we reconstructed seven bilateral ventral tegmental area (VTA)-subcortical pathways spanning direct mesolimbic projections and the hippocampus-VTA (HPC-VTA) loop regulatory circuit. Projection-endpoint analyses revealed long-axis and subfield organization of VTA-hippocampal pathways and superior-inferior VTA-nucleus accumbens (NAc) topography. Consistent with the HPC-VTA loop framework, neurite density covaried across mesolimbic pathways but not the arcuate fasciculus, demonstrating circuit-level organization. Nodewise analyses identified localized associations between pathway microstructure and delay discounting across VTA-hippocampal, VTA-NAc, and VTA-amygdala pathways. Together, these findings parse the human mesolimbic system into distinct pathways supporting adaptive behavior. The probabilistic atlas, MesoConnect, is freely available with tutorials to facilitate reconstruction and investigation of circuits in independent datasets.

The existence and mechanistic bases of sex effects on human subcortical anatomy remain unclear. To address these gaps, we combined structural neuroimaging of sex differences in subcortical volume and shape across three independent cohorts (combined N=2,886) with neuroimaging maps of sex chromosomal and endocrine effects from three clinical cohorts (combined N=313, including sex chromosome aneuploidies and isolated gonadotropin-releasing hormone deficiency). We found that human subcortical anatomy shows focal sex differences of moderate-strong effect size during young adulthood (male-biased surface area overlying lateral/ventral putamen and pallidum, nucleus accumbens, ventral posterior complex of thalamus, anterior hippocampus, and medial and central nuclei of amygdala; female-biased surface area overlying caudate head, posterior/dorsal pallidum, mediodorsal, lateral dorsal, lateral geniculate, and medial geniculate nuclei of thalamus, posterior hippocampus, and basolateral and lateral nuclei of amygdala). These regional sex differences show strong within-sample reproducibility at modest sample sizes (N[&ge;]200) and are replicable in a fully independent mid-adulthood cohort. Adult subcortical sex differences are already evident in nascent form by pre/early puberty (age 9 years), and the same holds for cortical sex differences. Comparison with neuroimaging data from rare disorder cohorts reveals a statistically significant spatial and directional congruence of subcortical sex differences with the anatomical effects of X-chromosome dosage in the striatum and effects of testicular hormones in the amygdala. Our findings provide a reference map for the fine-grained spatial patterning of sex differences in adult subcortical shape, demonstrating its replicability, establishment before pubertal progression, and spatial concordance with potentially causal sex chromosomal and gonadal effects.

Chinese is one of the living languages with the highest homophone densities, making homophone processing a critical component of language comprehension and orthographic decoding, as well as an ideal model for examining its cognitive mechanisms. This study assessed mono- and bi-character homophone selection function in 129 Mandarin- and Cantonese-speaking participants, including 66 cognitively healthy individuals, 17 with semantic variant (sv) primary progressive aphasia (PPA), 23 with the nonfluent/agrammatic variant (nfv) PPA, and 23 with the logopenic variant (lv) PPA. All PPA variants showed impairment in homophone selection, with most significant deficits in svPPA and lvPPA. ROC analyses demonstrated its strong diagnostic utility for distinguishing PPA from controls. Voxel-based morphometry revealed that homophone performance correlated with gray-matter volume in left anterior and inferior temporal regions, with variant-specific patterns observed. These findings described the potential role of homophone selection in diagnosing PPA syndrome in Chinese speakers, illuminate crosslinguistic influences on language impairment in PPA, and provide novel insights into the cognitive processing of homophones.

In aphasia, variability in lesion location and post-stroke network reorganization complicates efforts to identify the regions that support semantic versus phonological processing. In prior work, researchers have often examined these domains in separate patient cohorts, further limiting direct comparison of their neural substrates within aphasia. To address this gap, we investigated the degree of overlap and specialization within semantic and phonological networks in adults with and without aphasia and examined how neural activity relates to task performance and lesion profiles. Sixteen neurologically healthy adults and 15 participants with aphasia following left hemisphere stroke completed adaptive Semantic Matching and Rhyme Judgment tasks during functional near-infrared spectroscopy (fNIRS) recording. Changes in oxyhemoglobin and deoxyhemoglobin were measured within regions of interest, and activation was related to standardized task performance and damage to left ventral and dorsal language pathways in participants with aphasia. Controls showed largely distinct activation patterns consistent with prior fMRI findings: Semantic Matching preferentially recruited ventral stream regions, particularly left temporal cortex, whereas Rhyme Judgment more strongly recruited dorsal regions, including left inferior frontal and inferior parietal cortex. Participants with aphasia showed greater inter-individual variability and no clear group-level segregation of dorsal and ventral activity. Better Semantic Matching performance was associated with less left ventral stream damage and greater left dorsal activation, whereas better Rhyme Judgment performance was associated with lower right hemisphere ventral and dorsal activation. These findings suggest that post-stroke language performance reflects residual specialization, lesion-dependent constraints, and flexible recruitment of surviving tissue. Overall, fNIRS shows promise as an alternative to fMRI for language mapping in post-stroke aphasia, although replication in larger samples with more extensive cortical coverage is warranted.

Neuroimaging studies identify the superior temporal sulcus (STS) as an anatomical and functional target in autism spectrum disorder (ASD). Yet, the STS is a heterogeneous structure containing unique subregions, such as the superior temporal asymmetrical pit (STAP), which is a recently identified human-specific landmark. To date, no study has compared the STAP between ASD individuals and neurotypical controls (NTs) in a large sample. Here, we analyzed the mean intra-hemispheric sulcal depth and inter-hemispheric laterality of the STAP, as well as the anterior and posterior portions of the STS (aSTS, pSTS), in 100 NTs and 100 ASD male individuals (ages 5-18) using a discovery-replication design. Across samples, the STAP and pSTS showed rightward asymmetry, while the aSTS showed leftward asymmetry within groups. All features of the STAP, aSTS, or pSTS did not differ significantly between groups. However, pSTS laterality related to individual differences in social cognition, which was further supported by a meta-analysis across over 200 studies showing that neurocognitive functional differences associated with social cognition and ASD are restricted to portions of the STS outside of the STAP. Together, these findings provide fine-grained insight into prior studies showing coarse morphological differences of the STS between individuals with ASD and NTs, as well as open up new questions and new neuroanatomical targets of the heterogeneous STS in future studies of brain structure, brain function, and cognition in ASD, as well as other neurodevelopmental disorders targeting the superior temporal cortex.

Lesion Network Mapping (LNM) is a method used to map symptom-related circuits in the human brain. Recent work has raised concerns about the method and the specificity of LNM-derived networks, which often show high spatial overlap across unrelated conditions. Null models have been proposed to address the raised limitations1. Useful statistical tools in network neuroscience, they are reference procedures; in LNM, the normative connectome matrix C remains the source of the derived networks also under null modelling. Their broader role as a general solution for LNM therefore warrants careful evaluation. Here, we examine the proposed null models step-by-step and test them empirically across multiple datasets. We show that the proposed null-model correction effectively only replaces C with a transformed matrix C' (e.g., C' ~ C - mean(C)), one that is directly derived from C and retains much of its low-dimensional structure. Consequently, null-model-corrected LNM maps remain shaped by the same standard set of connectome-derived patterns and continue to show repetition and non-specificity across unrelated conditions. We discuss the scope and limitations of null models for LNM, showing that they can statistically benchmark LNMs, but that their application does not resolve the core limitation of non-specificity of LNM networks.

The posterior parietal cortex integrates somatosensory information with signals supporting action, spatial cognition, and multisensory processing, yet the fine-scale topographic relationship between the intraparietal sulcus (IPS) and primary somatosensory cortex (S1) remains poorly understood. Here, we used resting-state fMRI in 10 awake rhesus macaques to characterize functional connectivity between four S1 subregions (areas 3a, 3b, 1, and 2) and 223 single-voxel seeds densely sampling an IPS-centered parietal territory. Across S1 subregions and hemispheres, IPS-S1 connectivity exhibited a highly structured spatial organization characterized by systematic anteroposterior and mediolateral gradients and a prominent semi-oval topology following the geometry of the intraparietal sulcus. Connectivity was generally stronger toward the cortical banks and convexities and weaker near the sulcal fundus, a pattern that could not be explained by reduced local temporal signal-to-noise ratio. Unsupervised hierarchical clustering independently recovered this spatial organization and identified connectivity-defined subdivisions that only partially corresponded to classical cytoarchitectonic boundaries. Along the S1 mediolateral axis, prominent troughs in cluster-averaged connectivity profiles corresponded to several previously described somatotopic transitions, allowing to consistently map putative boundaries between tongue and face representations and between hand and trunk-leg representations. Based on these boundaries, individual IPS clusters exhibited connectivity with multiple S1 body-part sectors, with the pattern of connectivity across these sectors varying across IPS clusters and S1 subregions. Together, our findings reveal that IPS-S1 functional connectivity is organized by a large-scale sulcal topology while preserving fine-grained features of S1 somatotopic organization. These findings suggest that the functional organization of the IPS is shaped by multiple spatial dimensions, including large-scale sulcal topology and the organization of S1 subregions and body-part representations, which may together contribute to the functional specialization of individual IPS regions.

The glymphatic system is the waste clearance pathway in the central nervous system. Although its dysfunction has been associated with several neurological disorders, its role in multiple sclerosis (MS) remains poorly understood. This system is thought to be particularly active during sleep, yet the relationship between sleep metrics and MRI-derived markers of glymphatic activity has been scarcely explored. In this study, we aimed to perform a multimodal magnetic resonance imaging (MRI) assessment of glymphatic function in MS and investigate its association with sleep parameters, clinical disability, disease duration and MRI-derived structural markers. The study included 22 patients (34.3{+/-}10.1 years; 16 female) with relapsing-remitting MS and 26 healthy controls (32.2{+/-}5.7 years; 18 female), all of whom underwent 3.0 Tesla MRI. Glymphatic function was assessed using diffusion tensor image analysis along the perivascular space (DTI-ALPS) index and global blood-oxygen-level-dependent (BOLD) signal-cerebrospinal fluid (CSF) coupling. In addition, sleep micro- and macrostructural features, clinical disability, disease duration, along with mean diffusivity (MD), lateral ventricular, choroid plexus volumes and lesion load were evaluated. Patients showed reduced DTI-ALPS index, though this difference became only a trend after controlling for MD. No significant group differences were observed in global BOLD signal-CSF coupling. However, in patients, the weaker coupling was associated with longer total sleep time, lower arousal index, greater disability and lateral ventricle enlargement, whereas the DTI-ALPS index showed no significant correlations with any of the measures. Together, these findings challenge the assumption that DTI-ALPS can serve as a reliable marker of glymphatic function in demyelinating diseases such as MS, where disease-related alterations in white matter microstructure can confound diffusion-based measures. In contrast, gBOLD-CSF coupling showed associations with sleep, clinical disability and structural changes in patients, despite the absence of a significant group-level difference. This finding suggests that gBOLD-CSF coupling may capture individual variability in glymphatic-related dynamics relevant to sleep and the clinical expression of MS.

The superficial white matter (SWM), a thin layer of white matter beneath the cortex, contains a complex arrangement of short- and long-range fiber systems that support communication between cortical regions. These fiber populations include subcortical U-shaped fibers (U-fibers) that predominantly follow a tangential trajectory along the gray matter-white matter interface and long-range fibers that extend more radially through the deeper white matter to connect distant cortical regions. Here, we developed a surface-based framework to dissociate radial and tangential fiber components within the human SWM and investigated their microstructural properties and functional relevance using in-vivo high-field MRI. These components were quantified using apparent fiber density (AFD), fractional anisotropy (FA), and mean diffusivity (MD). We further contextualized SWM fiber architecture by examining depth-dependent microstructural profiles and their relationship with normative patterns of cortical organization and hierarchical specialization. The anatomical assessment of the proposed fiber orientation separation was supported by tractography-based bundle analyses, while the anatomical accuracy of tractography reconstructions was confirmed through comparison with post-mortem white matter dissection data. To assess their functional relevance, the fiber predominance map derived from radial and tangential AFD was related to spatial patterns of functional variability across the cortex. Regions with greater radially oriented fiber contribution exhibited more gradual changes in functional similarity across geodesic distance, whereas tangentially oriented fiber-enriched regions showed more rapid local functional differentiation. These findings reveal distinct structural and functional signatures of radial and tangential SWM fiber components and provide a framework for investigating their contributions to cortical organization and functional specialization. SIGNIFICANCE STATEMENTThe human cerebral cortex is organized into specialized regions that support diverse functions, but the structural principles underlying this organization remain incompletely understood. The superficial white matter (SWM), a thin layer of white matter beneath the cortex, contains intermingled short- and long-range fiber systems that support communication between cortical regions. However, these fiber populations have remained difficult to characterize because they are closely intertwined within the SWM. By developing a reproducible framework for mapping and quantifying distinct SWM fiber systems, this work provides a foundation for systematic investigation of this complex compartment. Integration of multimodal datasets, together with open sharing of analytical tools and standardized features, promotes transparency, reproducibility, enabling future studies of how superficial connectivity supports brain organization.

INTRODUCTION: The basal forebrain (BF) provides the principal cholinergic innervation to the cortex and the hippocampus and is among the earliest regions affected by neurodegenerative processes. Yet, whether its subnuclei differ in how their structural integrity changes with age and in their vulnerability to neurological and psychiatric disorders is unresolved. METHODS: Here, in 81,870 UK Biobank participants aged 45-87 years, we built age- and sex-specific normative models of grey-matter volume in the nucleus basalis of Meynert (Ch4) and the septal/diagonal-band (Ch1-3) compartments of the BF and derived individual deviation from these norms. We then tested the selectivity of their normative deviation profile across cross-sectional and longitudinal ageing, manifest neurological and psychiatric disorders, pre-diagnostic neurodegeneration, prospective disease incidence, and neuropsychiatric polygenic liability. RESULTS: The two compartments followed distinguishable age-related profiles across mid-to-late adulthood, which persisted when examined using longitudinal data. Deviation from norms was generally concentrated in Ch4 across neurological diseases and, in almost no condition, greater in the Ch1-3. The Ch4 deficit was further detectable years before diagnosis, was not substantially explained by common inherited genetic risk, and was linked with incident dementia and Parkinson's disease, whereas Ch1-3 carried no independent prospective structural vulnerability. DISCUSSION: These findings identify an anatomically specific age-related profile of the human BF, with structural vulnerability manifest and pre-diagnostic neurological disorders, which is not substantially explained by common inherited risk.

Bimanual coordination requires the control of each hand to account for movements of the other hand. To enable such coordination, it has been hypothesized that the brain has a representation of bimanual movements that supersedes a linear combination of the unimanual movements. To test for such a representation, we mapped the encoding of movement direction during unimanual and bimanual wrist movements in the human brain using functional magnetic resonance imaging (fMRI). For unimanual movements, we found that both the contralateral and ipsilateral movement directions were represented across most cortical motor regions. The representations were similar for mirror-symmetric movements across the hands, revealing an intrinsic (body-centric) code in premotor and parietal regions. For bimanual movements, the neural representations of movement directions could be explained by three components. First, the contralateral movement direction was represented as in the unimanual condition, whereas the representation of the ipsilateral movement direction disappeared. Second, in higher-order regions (rostral premotor and posterior parietal cortex) we found relatively elevated activity when the two movements were directed in unrelated directions. This activity was likely caused by the extra effort to deal with incongruent spatial targets across the hands. Third, and independent of this congruency effect, we found a non-linear interaction between contra- and ipsilateral movement that encoded the specific bimanual movement combination across sensorimotor, premotor, and parietal regions. This representation was suited to facilitate bimanual coordination and learning.

Structural differences observed in collision-sport athletes may reflect repetitive head-impact exposure, broader effects of athletic training, or pre-existing individual variation, making inclusion of a non-collision athlete group important for interpretation of any differences seen. Here, structural T1-weighted MRI data from 17 collegiate American football players, 16 cross-country runners, and 9 non-athletes were analyzed using FSL voxel-based morphometry (VBM), FIRST subcortical volumetry and vertex-wise shape analysis, and FAST whole-brain tissue segmentation. Spatial inference used permutation testing with threshold-free cluster enhancement and family-wise error correction, while scalar analyses used multiplicity-corrected procedures. FAST-derived absolute tissue volumes were also examined in a sensitivity analysis adjusting for FreeSurfer estimated total intracranial volume (eTIV). VBM showed recurrent brainstem-centered differences, with non-athletes showing greater local gray matter measures than both athlete groups, whereas football players showed greater brainstem measures than runners. FIRST vertex-wise analysis identified corrected differences involving the brainstem and selected subcortical structures, while no total FIRST structure volume survived false discovery rate correction. Unadjusted FAST analysis identified a CSF group effect, F(2, 24.38) = 11.46, p < .001, q = .001. After eTIV adjustment, the CSF group effect remained, F(2, 38) = 14.30, q < .001, and a white matter group effect was also identified, F(2, 38) = 9.59, q < .001; gray matter and total segmented tissue did not show corrected adjusted effects. These findings indicate that structural variation depends on anatomical scale, comparison group, and head-size adjustment and cannot be attributed uniquely to collision-sport exposure.